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ABSTRACT 

A self-consistent model for binary evolution devised by Gingold & Monaghan (1980) 
is used to show that two distinctly different types of behaviour are possible for close 
eccentric binaries. The model is based on a linear adiabatic normal mode analysis of the 
problem, which allows detailed examination of the transfer of energy from the orbit to 
the tides. We show that for most binaries, energy is exchanged quasi-periodically, with 
the system regulating itself so that the maximum tidal energy always remains small, and 
no circularization takes place. In contrast, for a range of eccentricities and periastron 
separations, chaotic behaviour prevails, with the eccentricity following a random walk and 
with the energy transferred to the tides during a single periastron passage being up to an 
order of magnitude larger than that transferred during the initial periastron passage. 

These results have important consequences for the study of tidal capture binaries. 
The standard model (see, for example, McMillan, McDermott & Taam 1987) assumes 
that the energy transferred to the tides during a periastron encounter is independent of 
the oscillatory state of the stars, and that the amount of tidal energy present at any time 
can be calculated using a formula which (accurately) gives the amount deposited after the 
first encounter (Press & Teukolsky 1977). The calculations presented here show that a 
self-consistent treatment is necessary in order to study the dynamical evolution of tidal 
capture binaries. We conclude that the tidal capture process would not be possible were 
it not for the existence of chaotic behaviour. 

1 Introduction 

The dynamical history of binaries formed by tidal capture (Fabian, Pringle, & Rees 1975) 
has been studied only for the first few orbits and for extremely close encounters (for ex- 
ample, Gingold & Monaghan 1980, hereafter GM, Kochanek 1992, Rasio & Shapiro 1991). 
The energy transferred during the capture process itself was first calculated accurately by 
Press & Teukolsky (1977) for ra = 3 polytropes and subsequently by Lee & Ostriker (1986) 
for n = 1.5 polytropes, these providing better models for the low mass main sequence 
stars found in globular clusters where tidal capture will be most effective. McMillan et al. 
(1987) considered the energy transfer in more realistic (non-polytropic) models. 

Several authors (for example, McMillan et al. 1987 and Ray, Kembhavi, & Antia 1987) 
have based their studies on the Press & Teukolsky (1977) calculation, and in particular, 
have extrapolated their results by assuming that the energy transferred during subsequent 
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periastron passages is independent of the pulsational state of the stars. The conclusion 
is then drawn that the orbit circularizes on a very short timescale (rcirc 10 yr). By 
then the oscillation energy is of the order of the binding energy of the polytrope, and 
since the timescale for the most energetic modes to thermalize is much longer than the 
circularization timescale (10'* — 10^ yr), the star necessarily responds by expanding, with an 
estimated order unity increase in the stars radius (McMillan et al. 1987). The conclusion 
is drawn that most tidal capture binaries result in collision or merger, given that the 
maximum periastron separation for the process to operate is about 3 stellar radii (for 
equal mass stars). 

The dynamical evolution of tidal capture binaries is clearly important in deciding their 
fate. 

The present work is based on an adiabatic linear normal mode analysis of the problem 
devised by GM. The self-consistent nature of the model (the total energy is conserved, 
allowing energy to flow freely between the polytrope and the orbit) allows us to follow 
the evolution of the system indefinitely (for as long as the energy is conserved to within 
some tolerance). While numerical methods such as Smoothed Particle Hydrodynamics 
(SPH - see GM) are good for examining extremely close encounters, particularly those 
near or at disruption, the present method allows us to examine relatively wide and hence 
stable orbits for the purposes of studying the formation and evolution of objects such as 
cataclysmic binaries, low-mass X-ray binaries, pulsar binaries, and black hole binaries, all 
of which (except, perhaps, black hole binaries) appear plentiful in the cores of globular 
clusters. 

The binary is modelled by a point mass and a polytrope of index n = 1.5, representing 
a compact object and a fully convective low mass Population II star respectively. The 
mass ratio is arbitrary. 

The equations of motion may be derived from a Lagrangian (GM) or directly by con- 
sidering the stars to be composed of iV particles and taking the continuous limit (Mardling 
1991). Either way affords a self-consistent treatment (and, of course, yields the same equa- 
tions), with the stars and orbit free to exchange energy in either direction. A normal mode 
analysis is then performed in which the velocity and density perturbations are assumed 
small, and a set of time-dependent ordinary differential equations is derived for the or- 
bit and the oscillation amplitudes of the modes. An arbitrary number of modes may be 
included. 

In Section 2, we review the Lagrangian model devised by GM. In Section 3, we cal- 
culate the orbital energy transferred to the tides after the first periastron passage of a 
parabolic capture orbit for a range of periastron separations and compare the results to 
those of Press & Teukolsky (1977). Section 4 considers the way the oscillation energy 
is distributed amongst the modes and the effects of resonance. In Sections 5 and 6, we 
present two distinct types of orbital behaviour; non-chaotic orbits which are quasi-periodic 
in eccentricity, and chaotic orbits for which the eccentricity follows a random walk. Sec- 
tion 7 considers for what range of initial periastron separations and eccentricities chaotic 
behaviour can be expected and plots a boundary between chaotic and non-chaotic orbits. 
Finally, Section 8 presents a discussion of the results. 

2 Equations of Motion 

The details of the derivation of the Lagrangian for the system may be found in GM. Here 
we present a summary. 

The motion of the system is referred to a non-rotating, non-inertial reference frame 
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with origin at the centre of mass of the polytrope. We assume that the polytrope is 
adiabatic and has no initial vorticity, and since the forces are conservative, the vorticity 
remains zero. The velocity of the fluid can therefore be written in terms of a potential, 
i.e. V = Vcj). If A is the position of the point mass relative to the centre of mass of the 
polytrope and the polytrope and point mass have masses Mi and M2 respectively, the 
Lagrangian L for the entire system is 
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Here p{r) is the fluid density at the point r, dr is a volume element and U = nnp^l^ is the 
thermal energy per unit mass, k being a constant and n the polytropic index. The first 
three terms respectively represent the orbital kinetic energy, the negative of the orbital 
potential energy and the negative of the interaction energy between the orbit and the 
fiuid, while the remaining terms represent the energy of the fiuid. Note that we have not 
included a term involving A-r, which GM later show to represent the linear momentum 
of the fiuid relative to its centre of mass, which of course vanishes. 

The Lagrangian may be scaled using the usual Chandrasekhar (1939) scaling: 



r=aX, <j) = {a^ I (3)^ , t = [it' , 



a- = {n+ l)KA(i-'^)/'^/47rG', GA/J^ = n/47r 



M = Aa^g andM2 = sMi, 
A is the central density so that -D(O) = 1 



p = \D, 
where 

and Q = 4tt J^" D{X)X'^dX, 

Xq being the scaled radius of the polytrope. As in GM, we retain the symbol A for the 
separation which will now be in units of a. We also drop the prime from the scaled time 
t' . Thus L = (Xa^ / (3)Lo, where the scaled Lagrangian is 
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The density is now written as _D = Dq + r], where Dq is the Emden density of an unper- 
turbed polytrope and r] as well as $ are assumed to be small. 

We now expand r] and $ in terms of the normal modes of vibration of the polytrope: 



k 
k 



(3) 
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where k = klni, 
aumit) = (-1)" 



Jkb 



i{t) = the amplitude of the density perturbation, 

^kl-m(^) amplitude of the velocity potential, r]ki and ([ki are the 

radial modes of oscillation as derived in GM (the relationship between these and the 
standard representation is presented in the Appendix), Yim are spherical harmonics as 
defined in Jackson (1975), (with Y/m^/'m"^^^ = ^ll'^mm') X, 9 and Lp are the spherical 
coordinates of a point in the polytrope and = YLV=i ^1^2 Y^m=-l- easily shown 
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that modes with / = 0, 1 are not excited. The k = 1 modes are the fundamental or /- 
modes, and the k > 2 modes are the p-modes. The g(-modes are not excited in an ra = 1.5 
polytrope (see, for example. Cox 1980, p238). 

The orthogonality condition for the radial modes is 

VkiA'i'X dX = ijJkihihk'^iv, (5) 



where Uki is the frequency of vibration of mode k. These are accurate to the same number 
of decimal places as those quoted in Lee & Ostriker (1986). To test the accuracy of the 
radial normal modes, we follow McMillan et al. (1987) by finding the second derivatives 
of (j)ki and riki numerically ((^^^ and 77^^ ) and comparing them with those obtained di- 
rectly from the differential equations. This is done by calculating a mass- weighted mean 
fractional error SkV- 
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Q Jo 

which is found to be between 10^^ and 10^^ for modes k = 1, ..,4 and / = 2, .., 10. 

Note that since the rotation (toroidal) modes (Unno, Osaki, & Ando 1989, pl88) are 
not excited, uj^i does not depend on the index m and so is (2/ + l)-degenerate. 

Expanding equation (2) to second order in $ and r], and using the orthogonality 
condition equation (5) (as well as that for the Y/m), the Lagrangian becomes 

Lo=2Y^s^ +^^A+''''^^ 27TT — A^^i — +2i^iKK-^kiKK)hu (7) 

k k 



where 



Tki= [""^kiX^^'dX (8) 



/o 

and Iki is defined in equation (5). The equation of motion of the orbit is thus 




while the amplitudes of oscillation are governed by 



with 



2, / nsQ Tu\ y4(vr/2,V>) 



«k = ^k/^fc/- (11) 



The second term on the right hand side of equation (9) represents the interaction between 
the tides and the orbit and is responsible for the deviation of the orbit from that of two 
point masses. The right hand side of equation (10) represents the forcing by the orbit of the 
oscillations of the polytrope. Equations (9) and (10) are integrated using a fourth-order 
Runge-Kutta scheme. 
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2.1 Initial Conditions 

We start the binary off at apastron (tp = —it, A = 0) or at a separation of about 40 stellar 
radii for capture orbits and supply the initial eccentricity, cq, and periastron separation, 
Ao, the binary would achieve as two point masses. For the first set of models, we start the 
polytrope off with zero oscillation energy (by^ = = 0)- We then consider models which 
start with a prescribed amount of oscillation energy. The initial values for the b]^ and b]^ 
are obtained from a run which is started with zero initial oscillation energy. 
In all models considered, we choose the mass ratio, s, to be unity. 

2.2 Energy and Angular Momentum 

The conserved (scaled) energy E and angular momentum / about the centre of mass of 
the system are shown in GM to be 

^ = 2YT1^ " A " ''''^ ^ 27TT A^^^ + 2 Z^(^k&k + ^kib],bi,)Iki, 

k k 

(12) 

and 

J = -^A^-Y,hK^mhi (13) 

k 

respectively, where i = \f—T. The first two terms on the right hand side of equation (12) 
represent the orbital energy, the next term represents the interaction energy between the 
orbit and the fiuid and the last two terms represent the oscillation energy of the fiuid. 
Note that since we take zero initial perturbation to the polytrope, the angular momen- 
tum remains in the direction perpendicular to the orbit. The first term of equation (13) 
represents the orbital angular momentum, while the second term represents the angular 
momentum transferred to the tides, which remains small relative to the total angular 
momentum. 

The total energy and angular momentum are conserved to at least one part in 10® in 
all models considered in this paper. 



3 Comparison with the Press and Teukolsky Model 

Press & Teukolsky (1977, hereafter PT) calculated the amount of energy deposited in a star 
after the first periastron passage of a capture orbit. We first show how their calculation 
is related to ours, and then compare our model to theirs by repeating their calculation. 

In the present units, the rate at which energy is transferred from the orbit to the 
extended star is given by 

f^|z)v.VJ.dX. ,14) 

where v is the velocity of a fiuid element and U is the gravitational potential of the point 
mass, given by 

Assuming that the orbit remains parabolic, the amount of energy deposited after one pe- 
riastron passage is then J^^{dE / dt)dt. PT calculate this quantity by writing the velocity 
in terms of the normal modes of the star, taking the density to be the unperturbed density 
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and assuming that the orbit remains parabolic. In terms of the present model, one may 
show that in fact, 

^ = ^E(W + ^'«)W2. (16) 

which is the rate at which at which oscillation energy is deposited or removed from the 
polytrope. If we remove the orbit-oscillation interaction term from equation (9) so that 
a capture orbit remains parabolic, we find the relative difference between our model and 
that of PT is 0.006% for a periastron separation of 2 stellar radii, (polytropic index 1.5) 
increasing to about 0.1% for a periastron separation of 3 stellar radii. On the other hand, if 
the interaction term is left in, our model predicts about 4% more energy will be deposited 
after a 2 stellar radii encounter, and this drops to 0.7% at 3 stellar radii. 

The latter differences between the predictions of the two models are insignificant for 
the first periastron passage, but we show in the following section that the PT model can 
underestimate the energy transfer during a single subsequent passage by up to an order 
of magnitude if the resulting orbit turns out to be chaotic. On the other hand, we show 
that for non-chaotic orbits, the energy exchange between the orbit and the polytrope is 
quasi-periodic, so that energy is not steadily added to the polytrope at each periastron 
encounter, as was previously thought. 



4 Mode Energy Distribution and Resonance Effects 

Although in this type of study a few p-modes are usually included in calculations, they 
in fact carry little energy. For example, after the first periastron passage of an initially 
parabolic orbit, the pi modes (k = 2) carry less than 0.04% of the total oscillation energy 
when the periastron separation is 2Xo, and only 0.0001% when it is SXq. This remains 
true after subsequent periastron passages. 

For the /-modes (k = 1), after the first periastron passage of an initially parabolic 
orbit with Aq = 2Xo, the energy residing in the / = 3 and / = 4 modes is about 13% and 
2% respectively of that residing in the / = 2 mode, while for such an orbit with Aq = 3, 
the corresponding figures are 5% and 0.5%. This does not necessarily remain true for 
subsequent orbits. If the orbit happens to be chaotic (see below), it can happen that the 
/ = 3 mode carries more energy than the / = 2 mode, and this can become important when 
a capture orbit (cq = 1) 'tries' to avoid ionization when in the presence of other stars (see 
Mardling 1994, hereafter Paper II). On the other hand, if the orbit is not chaotic, the 
above distribution of energy persists. 

Finally, as far as the m-modes are concerned, most of the energy in the /-modes resides 
in the / = m modes. After the first periastron passage of an initially parabolic orbit with 
Aq = 2Xo (Aq = 3Xo), the energy residing in the / = m modes is 97% (99.5%) of the 
total oscillation energy, although the distribution is not as severe for less eccentric orbits. 
As we will show in a future publication (Mardling 1994b), when non-linear mode-mode 
interactions are taken into account, it is possible for high /-modes to be excited, and for 
these /-modes it is the m = modes that are most excited (compared to other m-modes 
with the same value for /). 

The tides of stars in close binary orbits are principally raised by the coupling of the 
orbit and the / = m modes. Consider the term which forces the oscillations in equation 
(10). Recall that y/„(7r/2, t/') = Q„P;™(0)e™'^ (Jackson 1975), where cim is a constant 
and P™ is a Legendre function of degree / and order m. Clearly, resonance effects occur 
if it happens that nitp ^ uj^i. This is in fact true for / = m when an orbit is at periastron, 
as can be seen in Figure 1 which plots curves in (Ao,eo) space along which Itpp/ujki = r 
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Figure 1: Resonance curves in tlie (Aq, eo) plane, along wliicli Itp/ujki = r for r = 1 and 
r = 1/2 and / = 2, 3, . . . , 6. Here and in tlie figures which follow, Aq is measured in units 
of the scaled stellar radius, Xq- 

for various /, and for r = 1 and r = ^ {ipp may be written in terms of Aq and cq by using 
the equation for an unperturbed orbit). Here and elsewhere, the subscript p refers to the 
value of the orbital variable at periastron. 

5 Non-Chaotic Orbits 

In order to examine the dynamical behaviour of these systems, we plot the orbital ec- 
centricity against the number of periastron passages. The eccentricity actually varies 
continuously, with maximum changes occuring near periastron. Nonetheless, the approx- 
imate eccentricity can be obtained via e = (1 — r)/(l + r), where r is the ratio of the 
apastron separation to the following periastron separation. Figure 2 shows the evolution 
of the orbit with initial conditions (Aq = 3.2Xo, eo = 0.7). The / = 2 modes are included, 
and the solution is not very sensitive to the inclusion of more modes: the deviation from 
this model by one which includes modes up to / = 4 is barely discernable after 20 orbits. 

Energy is periodically extracted from the orbit and then replaced. The direction of flow 
of energy depends on the relative phase of the orbit and the oscillations at periastron, with 
the long term behaviour exhibiting a beating effect. This is in contrast to the assumption 
made by McMillan et al. (1987) and others that energy is extracted from the orbit at each 
periastron encounter so that the orbit steadily circularizes. 

Figure 3 shows the advance of the apsidal line for this model, which is 26% greater than 
the classical prediction (see, for example, Schwarzschild 1965). This difference reduces to 
zero as the separation increases, as indicated in Figure 4 for an eccentricity of 0.4 (all 
orbits calculated for this figure are non-chaotic). 

The example in Figure 2 is representative of the general behaviour of most non-chaotic 
orbits. Other systems, such as that shown in Figure 5 for the case (Aq = 2.87Xo, eo = 0.2) 
exhibit a longer period and a larger amplitude, indicating that more energy is extracted 



8 



0.7002 




0.6992 



100 150 
periastron passage 



200 



Figure 2: A non-chaotic orbit for wliicli tlie eccentricity versus periastron passage ex- 
liibits beating. Note tliat tlie total cliange in eccentricity is small. Here, Aq = 3.2Xo 
and eo = 0.7. 
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Figure 3: The advance of the apsidal line for the model shown in Figure 2. The rate of 
advance is constant for non-chaotic orbits. 
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Figure 4: Comparison of tlie classical and dynamical calculations of the apsidal advance 
as a function of initial periastron separation, Aq (for non-chaotic orbits only). The 
fractional difference decreases with increasing Aq. Here, e = 0.4. 
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Figure 5: A non-chaotic binary with a long period in the eccentricity. Aq = 2.87Xo and 
eo = 0.2. 
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Figure 6: Two initially close chaotic orbits for which (Aq = 2.8Xo, cq = 0.8) and 
(Ao = 2.8OOOIX0, eo = 0.8). The orbits diverge after the 4"^ orbit. The total change in 
eccentricity is not bounded as it is for a non-chaotic orbit (but see Paper II). 



from the orbit than in the typical case, but the important features to keep in mind are that 
the average eccentricity is constant and that the amplitude of the variation is of the order 
of the change in eccentricity after one periastron passage for most cases, so that the tidal 
energy always remains small and no circularization takes place (at least in the absence of 
dissipation). These last two points have important consequences for the evolution of tidal 
capture orbits. This aspect will be examined in Paper II. 

In Section 7 we will consider for which parameter values Aq and cq this type of be- 
haviour is found, when we plot the boundary between chaotic orbits and non-chaotic 
orbits. 

6 Chaotic Orbits 

Figure 6 compares the evolution of the two orbits with initial conditions (Aq = 2. 8X0, eo = 
0.8) and (Aq = 2.80001Xo, eo = 0.8), while Figure 7 compares the apsidal advance. The 
behaviour is in stark contrast to that of non-chaotic orbits and clearly depends sensitively 
on the initial conditions. It is also extremely sensitive to other changes such as the inclusion 
of more modes, and changes in the timestep used to integrate the equations (and even the 
machine one performs the calculations on!) 

The eccentricity varies over a wide range of values, and any one change in eccentricity 
after a periastron encounter can be up to an order of magnitude more than the initial 
change. 

Other ways of establishing whether or not a solution is chaotic include determining 
the Lyapunov exponents of the system (see, for example, Rasband 1990), and plotting a 
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Figure 7: The apsidal advance for the model shown in Figure 5. The rate of advance is 
not constant for chaotic orbits, and can even be negative. 



Poincare surface of section (Rasband 1990). The number of degrees of freedom is large 
for this system, even when only the / = 2 mode is included. The author has developed a 
method for calculating the Fyapunov exponents of large systems efficiently (Mardling 1991, 
1994c), based on a method devised by Shimada & Nagashima (1979), and independently 
by Benettin et al. (1980). Applied to the present system with just the / = 2 modes included 
(see Mardling 1991), we find the largest Fyapunov exponent for a typical chaotic orbit to 
be about 0.01. 

6.1 Surface of Section 

A traditional way of visualizing a system's transition to chaos is to plot its Poincare surface 
of section. Normally this is only possible when the system has 3 degrees of freedom, which 
occurs for 3 dimensional systems and higher dimensional systems for which the number of 
integrals of the system reduces the degrees of freedom to 3. 

If we include only the / = 2 /-modes, it is possible to plot a surface of section for the 
present system, even though it has 6 degrees of freedom (after the energy and angular 
momentum integrals are taken into account). 

To see how we may do this, let us first recall how a Poincare surface of section is 
drawn for a system with three degrees of freedom. For example, for the four dynamical 
variables, x, y, x and i), with integral f{x,y,x,y) = E, where E is a constant, we might 
plot X versus x every time the solution curve crosses the y = plane, with i) being given 
by the integral. We would then observe whether successive points tended to lie on curves, 
indicating the existence of a further integral, or even whether the points were mapped to 
a finite number of points, indicating that the system was completely integrable and that 
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the solution was periodic. Chaotic motion would be indicated if the iterates tended to be 
space filling. 

This method does not normally work for more than three degrees of freedom because 
in order to examine the behaviour of the orbit as it crosses a two dimensional subspace, 
we require, for example, several variables to pass through zero coincident ally which is 
unlikely to occur in most situations. As it happens, a judicious change of variables makes 
this possible for the present model when we only include the / = 2 /-mode. 

We imagine that our set of dynamical variables is replaced as follows^ : 

{V', A, V^, A, (5k, b*^), (5k, ^k)} ^ e. A, Xk, ^kl, (17) 

where E is the total energy given by equation (12), / is the total angular momentum given 
by (13), e is the orbital eccentricity, Xk interaction energy between the orbit and 

the k*^ mode and _Ek oscillation energy in the k*^ mode, with k = (1,2,0), (1,2,-2) 

and (1,2,2). Xk -^k defined as follows: 

Ik = -nsQ^^^iby,Y,r, + 5^y4), (18) 

i^k = (W + ^'/W)4//2. (19) 

At apastron, A = 0, and for sufficiently high eccentricity, Xk — that the orbit 
and oscillations are essentially decoupled. Thus e and the E]^ are approximately constant. 
We also find that for non-chaotic orbits, the energy in the m = modes drops to zero at 
apastron, so that since E and / are constant, we are left to plot the energy in the (1,2,2) 
mode against the eccentricity. 

Figure 8 shows a series of surface of section plots for models with cq = 0.6 and Aq = 
2.845Xo to Aq = 2.848Xo. Figure 9 shows the final plot in Figure 8 for later times, while 
Figure 10 plots the eccentricity against periastron passage for this model, showing how 
the orbit passes through two non-chaotic phases, seen clearly in Figure 9. 



7 The Chaos Boundary 

We can do a systematic search for chaotic systems in the (Aq, cq) parameter space by 
comparing initially very close orbits and seeing if they diverge after the first few orbits. In 
this way, we can dertermine the 'chaos boundary' above which chaotic behaviour ceases. 
Figure 11 shows the results of this search. The sections of the curves corresponding to 
small eccentricities are not of much practical value: models starting with these values 
for Aq and cq disrupt when non-linear terms are included in the equations of motion. 
Intermittant chaotic behaviour is observed above the curves shown: there must exist a 
curve beyond which all chaotic behaviour has ceased. Nonetheless, the curve shown has 
practical value, as we show in Paper II, where we consider capture orbits. 

All models calculated so far have been started with zero initial oscillation energy. We 
can also plot chaos boundaries for models with specific initial oscillation energies. These 
also prove to be useful for examining capture orbits. These are shown in Figure 11 as the 
terminating curves: they terminate at the point in (Ao,eo) space which corresponds to 
zero total energy. Thus, of course, the 'zero oscillation energy' chaos boundary terminates 
at a point on the line e = 1. 



^Recall that ^k ^^'^ ^k independent. 
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Figure 8: A route to chaos: the Poincare surfaces of section for models with eg = 0.6 
and Ao between 2.845 and 2.848. Read from left to right, top to bottom. 
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Figure 9: Chaotic and non-chaotic motion intertwined: this shows the final plot in 
Figure 7 (the lower right-hand corner of this figure) for longer times. 
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Figure 10: Chaotic and non-chaotic motion intertwined: the eccentricity versus perias- 
tron passage for the model shown in Figure 8. 




Figure 11: Chaos boundaries. A chaos boundary labeled C'x is calculated by starting 
the polytrope off with an amount x of oscillation energy. These curves must terminate 
at a point corresponding to zero total energy. 
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8 Discussion 

We have shown that two distinctly different types of behaviour are possible when a fully 
self-consistent treatment of the dynamical evolution of close eccentric binaries is con- 
sidered. For non-chaotic orbits, the eccentricity of the orbit is quasi-periodic, with the 
maximum tidal energy always remaining small. In contrast, the amount of energy trans- 
ferred to the tides in a chaotic orbit can be substantial. The process should not depend 
on the structure of the star being polytropic, rather it is the mass distribution and hence 
the ability of a star to store tidal energy which is important. 

One fascinating aspect of this work is that the range of periastron separations for which 
the tidal capture process is possible (which depends on the local velocity dispersion), coin- 
cides almost exactly with the range of periastron separations for which chaotic behaviour 
occurs. This is true for stars with mass distributions something like an ra = 1.5 polytrope. 
For more centrally condensed stars, the range of periastron separations which make tidal 
capture possible reduces (McMillan, Taam & McDermott 1990), while the chaos boundary 
can be expected to drop (the tides cannot store as much energy). It remains to be seen 
whether these ranges coincide as they do for the present case. 

In conclusion, we may say that the tidal capture process would not be possible were it 
not for the existence of chaotic behaviour. 
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APPENDIX 



Relationship Between the Standard and Gingold & Mon- 
aghan analyses 

The standard normal mode analysis (Unno et al. 1989, p90) involves a displacement vector 
^ which can be expressed as a sum of radial and poloidal components: 

^ = J2 '^os(ukit) {^^i(r)er + ef/(r)rv) Y^O, ^) (20) 
k 

which is scaled in natural units (but see later). The velocity potential used in the present 
analysis may therefore be related to ^ via the (unsealed) velocity of the fluid: 

where M* and i?* are the mass and radius of the star respectively. The velocity displace- 
ment may be written in terms of normal modes as (GM) 

$ = - ^ Mujkit)MX)YUO, ^) (22) 
k 

so that 

= - Y,^m{ojkit){^ex + ^Xv\YUe,^). (23) 
k 

Thus from equations (20), (21) and (23), we have 

,R_ /47rXoy/2 1 d<i,kl .... 



and 



\ nQ J LOki X 



Note that here ^ is scaled such that (Unno et al. 1989) 

l""' KO^k-^k'dr = ^kk' l""' P(r) [{0 + 1(1 + 1) {fkl) 



r'^dr, (26) 



which is used by MMT, while it is more common (PT, Lee & Ostriker 1986, Ray et al. 
1987) to normalize ^ such that 



rV(0^k-^k'dr = ^kk'- (27) 
Jo 

This can lead to some confusion. Note also that the factor p{r) in the integrand of equation 
(27) is missing from this expression in Lee & Ostriker (1986). 

Since ^'^ is related by an algebraic expression to the pressure (and hence density) and 
gravitational potential perturbations, we may relate this to 77 in a simple way (see Unno 
et al. 1989 and GM). 
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